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ABSTRACT 


Loads methods currently being used to design Shuttle Orbiter payloads are 
summarized. Methods used for the design of payloads launched by expendable 
launch vehicles are described in historical perspective. Experiences gained 
from expendable launch vehicle payloads are used to develop methodologies for 
the Space Shuttle Orbiter payloads. The objectives for the development of a 
new methodology for the Shuttle payloads are to reduce the cost and schedule 
for the payload load analysis by decoupling the payload analysis from the 
launch vehicle to the maximum extent possible. Methods are described for 
payload member load estimation or obtaining upper bounds for dynamic loads, as 
well as load prediction or calculating actual transient member load time, 
histories. 
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I. INTRODUCTION 

The advent of the space shuttle as the major transportation system for 
U.s. space operations in the 1980's and beyond has renewed the interest in 
developing improved methods for payload dynamic load prediction. 

Consistent with the objective of the space shuttle system to achieve a low 
cost launch capability, cost-effective methods for obtaining payload 
design loads are being sought. To achieve these objectives, the Dynamic, 
Acoustic, and Thermal Environments (DATE) program has been formed by the 
NASA Office of Aeronautics and Space Technology (OAST). The DATE working 
group consists of representatives from NASA centers interested in 
structures technology and is chaired by the Goddard Space Flight Center 
(GSFC) . The goals of DATE are the development and validation of improved 
technology for the prediction of dynamic, acoustic, and thermal 
environments and the associated payload responses in cargo areas of large 
reusable vehicles. To achieve this, extensive research-oriented 
environmental measurements on early shuttle flights are to be made. These 
data will then serve as the basis for developing better load prediction 
methods which are to be confirmed and refined as additional flight data 
become available. Ultimately these data will lead to improved spacecraft 
structural design criteria and test criteria [l],r2]. — 

This report summarizes the loads methods for the Shuttle Orbiter payloads 
which are now being proposed for future usage by the payload community. 

A historic background of payload design methodology for expendable 
boosters is included for completeness. While this sugary emphasizes the 
experiences at the Jet Propulsion Laboratory (OPL), it is representative 
of aerospace practice over the last 20 years. In the context of loads 
analysis, the term "spacecraft" is frequently used in this report to mean 
the payload structure carried by the launch vehicle or expendable booster 

Extensive experience has been gained in the past at OPL in the structural 

design of such spacecraft as Ranger, Mariner, Surveyor, Viking, and 

Voyager. All these were launched using expendable boosters. 
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The Galileo spacecraft currently being developed at JPL is the first 
interplanetary spacecraft to be launched using the Space Transportation 
System (STS). Over the years spacecraft structures have been designed to 
various criteria. In the early days of the space program the structure 
was typically designed to survive a sine vibration test at input levels 
prescribed by the environmental requirements engineers. Later as 
computers and finite element methods became prevalent, structures were 
designed to survive an analytically predicted flight load. Lately methods 
were developed at JPL to calculate the expected upper bounds of the member 
loads. 

Since it is anticipated that ultimately the STS is to be used much like a 
cargo airplane, such that only minimal or no system preflight loads 
analysis will be required for high reliability missions, it is important 
to develop several tools for loads predictions. The methodology to be 
used for a particular payload will be determined by the size, complexity, 
and cost of the payload, and thus by the degree of dynamic interaction 
with the STS. The degree of dynamic interaction between the payload and 
the launch vehicle is much more important for a reusable launch vehicle 
due to the emergency landing loads, a condition wherein the 
spacecraf t-to- 1 aunch vehicle weight ratio is higher than for typical 
launch events of expendable boosters. 

A classification of STS payloads as to cost, size, and criticality as 
related to the detail of predicted analysis levels is given in Table 1-1. 

The purpose of this report is to summarize the STS loads methodology 
currently being used by the payload community and the methods under 
development to simplify the STS loads prediction. Wherever possible, the 
proposed methodology will be evaluated using actual design loads obtained 
for the Galileo (GLL) and the International Solar Polar Mission (ISPM) 
spacecrafts. 

This summary is intended as a first step in the development of loads 
methodology as part of DATE, and will be updated as more experience is 


- 3 - 


Table 1-1. Classification of STS Payloads 



Category I 

Category II 

Category III 

Cost 

Low cost, 
detailed 
transient anal- 
ysis too severe 
a cost penalty. 

Moderate cost, some 
loads analysis 
justifiable . 

Highly expensive, large 
cost penalty for 
structural failure , 
detailed loads analysis 
justified . 

Size 

Small/light, no 
significant 
effect on 
orbiter, can 
afford wt. 
penalty. 

Some impedance 
effects at orbiter 
interface in low- 
frequency range. 

Large payload/orbiter 
dynamic interactions, 
may effect orbiter/ 
external tank dynamics 
at launch. 

Criticality 

No critical 
launch window, 
can be 
retrieved/ 
repaired, 
numerous 
launches . 

Some duplication of 
payloads, con- 
straints on launch 
windows, minimal 
repairability 

Critical launch win- 
dows, no retrieve/ 
repair capability, 
single backup 
spacecraft. 



EXAMPLES 



Communication 
satellites , 
small science 
payloads, earth 
resources, wea- 
ther satellites. 

High energy/ 
synchronous, DOD 
payloads, manned 
orbital lab. , large 
space telescope. 

Planetary missions, 
comet rendezvous, 
out -of -ecliptic . 






- 4 - 


gained in applying these methods and comparing the results to flight 
data. This report is mainly concerned with the method of obtaining loads; 
however, the modelling techniques, analysts/test correlation of the 
mathematical and test philosophy will be discussed to the extent that they 

relate to loads methods. 
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II. EARLY LOADS METHODOLOGY FOR EXPENDABLE BOOSTERS 

The structural analysis and design engineer has been searching for a 
systematic method to obtain dynamic loads data since the design of 
Explorer I in 1957, long before spacecraft structural design criteria and 
spacecraft environmental requirements were formalized. At that time, 
spacecraft designers lacked the experience and technology of the aircraft 
industry both in the area of loads prediction and dynamic testing. Thus, 
the first Explorer was conservatively designed by a combination of 
engineering judgement and requirements to de-couple the spin-frequency 
from structural frequencies. Structural modeling for the Explorer 
satellite was limited to simple back-of-the-envelope type of 
calculations. The vibration test requirements were derived from existing 
military specifications for the Sergeant and Corporal Missiles. 

As the lunar and interplanetary spacecraft became larger and more complex, 
the design for a dynamic environment became more important. 

A. The Ranger Lunar Spacecraft 

Ranger was the first JPL project for which an attempt was made to 
design a spacecraft for a dynamic environment. The Ranger spacecraft 
weighed approximately 340 kg (750 lbs) and were launched__using the 
Atlas/Agena Space Vehicle. 

1. Structural Design Criterion _ _ 

Although the fundamentals of the dynamics of space vehicle 
systems were understood at the commencement of the Ranger 
project, they were not documented until much later 13). Since 
there was a lack of flight data, the structural design 
philosophy for Ranger was basically governed by an environmental 
specification for a system vibration test. An implied design 
goal was a minimum frequency criterion for . the cantilevered 
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spacecraft which was negotiated with the launch vehicle 
contractor. 


2 . Dynamic Analyses 

Since the computer tools were limited at the time Ranger was 
designed, only a limited amount of static and dynamic analyses 
were performed. The modeling was limited to a finite element 
representation using truss and beam elements C4J. The dynamic 
analyses consisted of estimating load factors for various masses 
and then solving an equivalent static problem. Thus the dynamic 
ana ysis relied heavily on engineering judgement, and was mainly 
aimed at successfully passing the vibration test. 


3. Structural Testi 


ng 


Structural testing on Ranger consisted of developmental testing 
and qualification testing. 

a. Developmental Testing 

Extensive modal testing using multiple shakers was 
performed on the Ranger spacecraft. The data obtained in 
these tests consisted of frequencies, mode shapes, modal 
damping and a limited number of force coefficients. 

Although data obtained in these early modal tests were not 
used to correlate with analysis, these data proved very 
valuable in re-evaluating design load factors and assessing 
the dynamic responses expected during qualification 
testing. On later Ranger spacecraft, the modal data was 
used to redesign the spacecraft to avoid coupling with 
space vehicle modes. The Ranger modal data proved very 

useful for the Surveyor spacecraft as will be discussed 
later. 








h* Qualification Tost, inn 

Iho qualification test requirements for the early Hanger 
spacecraft prescribed a three axis system vibration test to 
a level defined by the project environmental requirements. 
The input frequency spectrum was flat. This resulted in 
over testing and required some redesign at the cost of 
structural weight. Later in the project, a limitation on 
the vibration input levels based on the structural 
capability of the Aqona adapter was implemented. The 
justification was not. to test the spacecraft to levels 
higher than what the support structure could react, thus 
ensuring that in flight t he launch vehicle support 
structure would fail before the spacecraft would. 

4. Flight Instrumentation 

Flight instrumentation on the Ranger program was minimal ;„d was 
mainly aimed at defining a torsional pulse which had bean 
observed in flights of Atlas vehicles at the time of Booster 
Engine Cutoff (RECO). Two accelerometers were flown on later 
Ranger spacecraft to define this pulse. These flight 
measurements did little to affect the design of the Ranger 
spacecraft other than to provide frequency information, but 
pt ovod to be very valuable for the Surveyor program. 

fl. The Surveyor Lunar Spacecraft 


Surveyor was the first JPL project in which a limited attempt was 
made to use loads analysis methods to either design the structure 
th], Lb ) or define the qualification testing [7j. The Surveyor 
Spacecraft weighed approximately ?80 kg (680 lbs) and was launched 
using the Atlas/Centaur Space Vehicle. 


- 8 - 


1. Structural Desiqn Criteria 

The Surveyor structural design requirements were dominated by 
the expected 3 o lunar landing loads. The launch loads were 
generally expected to be lower. For the launch condition the 
structure was basically designed to survive a sinusoidal 
vibration test, the input level of which was intended to 
envelope the expected launch acceleration. A frequency 
criterion was not imposed on the spacecraft design. 

2. Dynamic Analyses 

Analytical tools for structural design were becoming readily 
available at the time the Surveyor analysif was performed. Not 
only were the finite element capabilities increasing l3l , but 
also the computer core availability allowed the solution of 
larger structural problems. 

Dynamic landing simulations were performed tor the purpose of 
assessinq landing stability [9], CIO], [11] and landing loads 
[12]. The landing analyses assumed a rigid spacecraft body 
using an elastic simulation of the non-linear landing 
structure. The equations of motion were integrated using a 
fourth order Runge-Kutta procedure. To obtain spacecraft 
landing loads, the interface forces between the landing struts 
and the spacecraft calculated in a landing simulation were 
applied to the elastic spacecraft, simulated by experimentally 
determined normal modes. This procedure was not entirely 
consistent because the stability simulation program did not 
account for elastic responses of the spacecraft. However, in 
this simulation an approximation with deviations generally in 
the conservative direction was expected. Several redesigns, 
especially in the antenna solar panel positioner, were performed 
based on the results of the landing load investigation. 


. . iw mmmmfmm 


II W9 m WWWWP 




- 9 - 


for the simulation of launch I nails, a I i mi toil amount, of analysis 
was performed using forciiui functions ostimatod from flight data 
measured during the Atlas/Agona/Ranger boost phase. These 
analyses were limited to estimates of spacecraft response to a 
torsional forcing function since that was the only flight data 
available [6], Similar types of analyses using the same flight 

data were later performed for the GSFC Orbiting Geophysical : 

Observatory (OGO-E) Spacecraft [ 13] and the Orbiting 
Astronomical Observatory (0A0 A-?) Spacecraft [14]. 

Structural Testinq 

Two types of structural testing for the Surveyor Spacecraft were 
performed: developmental testing and qualification testing. 

a. Developmental Testing 

Extensive modal testing with multiple shakers was performed i 

both by the Hughes Aircraft Company (MAC) and by JPL. The 

data obtained in these tests were used both for modal 

correlation and for obtaining dynamic data used directly in 

the landing loads analysis L 1?]. The modal test results 

for the boost configuration were used directly in the 

torsion launch load analysis! Tv], [.fil Structural transfer 

functions were measured for the purpose of assessing flight 

stability during the terminal descent phase [is],. 

b . Qu a 1 if ica t. ion Test inn 

Qualification testing was conducted for the landing loads 

and the launch loads. For the landing loads a prototype i 

spacecraft with all components in flight-like configuration 

and operation, including the telemetry link, was ■ 

type-approval tested for landing in three drop tests. This j 

was done in order to establish confidence in structural and 
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functional survival of the spacecraft when subjected to the 
dynamic landing environment. 

For the launch environment a sinusoidal vibration test on 
the system level was performed, allowing the input levels 
to be suppressed or "notched" at frequencies where the 
response of certain components such as the solar panel 
positioner exceeded the design loads as determined from a 
landing loads analysis. 

The lack of simulation of the boundary conditions during a 
sinusoidal vibration test was recognized during both the 
Ranger and Surveyor Spacecraft [16]. It was shown 
analytically that severe over- or undertesting can occur. 

A rocking test table allowing the proper truncation and 
rotation inputs was proposed. This mode of testing was, 
however, never implemented due to cost and complexity. 

4. Flight Instrumentation 

Flight instrumentation on the Surveyor Spacecraft was limited to 
strain gauges on the landing gear to measure landing forces. 

Data obtained from these measurements compared favorably to 
analytical predictions [11]. There was no flight 

instrumentation flown to measure spacecraft responses during the 
boost phase. 

C. The Mariner Interplanetary Spacecraft 

The Mariner Spacecraft series was the first JPL project for which 
m-house loads analyses were performed for some of those flight 
events for which flight data was available. The Mariner spacecraft 
ranged in weight from 245 kg (540 lbs) to. 980 kg (2150 lbs) and were 
launched atop the Atlas/Centaur Space Vehicle. 
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1 . 


Structural Design Criteri- 


on 


2 . 
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A(wi = H(<o) F(m) 


(II-l) 


where A(w) is the Fourier transform of the response acceleration 

F(u>) is the Fourier transform of the forcing function 

H(co ) is the complex structural transfer function for the 
space vehicle, and 


on's the circular frequency, rad/sec. 


It is easily seen that if A(to) is obtained from a previous 
flight and the transfer functions for the previous and the new 
space vehicle are Hq(w) and H^(w), respectively, the new 
spacecraft response A(w) can be calculated from 

a n M = H n (u)) h 0 h a A q M (II-2) 

assuming that F (to) is invariant. 

The above method was used to estimate Mariner loads from 
previously flown spacecraft, both Mariner and others. It should 
be noted that while the method is mathematically exact there are 
several practical limitations. The accuracy of this method 
depends on obtaining phase coherent data with a minimum of 
telemetry noise. Furthermore, the method is inherently 
sensitive to the fidelity of the structural model. A difference 
of computed modal properties with the observed frequency 
response will cause a spurious response in the inverse process 
of Eq. ( 1 1-2) . 

The frequency domain approach was used in the Mariner Program 
both on the digital and analog computer. In addition to flight 
load estimation, this technique was also used for determining 
the shaker armature overturning moment during the low frequency 
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longitudinal 

requirements 


vibration test and for determining the control 
during the forced vibration test. 


3. Structural Testing 


On Mariner, several unique development tests were conducted in 
addition to the qualification vibration test. 

a. Developmental Testing 

Frequencies, mode shapes and damping characteristics for 
the cantilevered Mariner Spacecraft were determined using 
mmiple shaker sine dwell tests. The modal data were used 
analytical model correlation. These test data were 
also used directly in the dynamic analyses. The modal 
survey showed that the propellant played a very important 
role in the proper representation of the dynamic model. 
Hence special dynamic testing using a single tank was 
performed CM. The objective of this testing was to study 
the dynamic behavior of the fluid, specifically the 
effective weight and the damping. 

b. Qualification Testing 

The major structural elements of the Mariner Spacecraft 
were qualified by a static test using strain gauged 
members. The static loading conditions were determined 
from a combination of predicted flight loads and expected 
oads during the qualification vibration tests. Forced 
vibration testing was performed on the Mariner Spacecraft 
system to levels specified by the environmental requirement 
group. The input levels were suppressed in order not to 
exceed design levels for primary structure. 
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4. Flight Instrumentation 

Several of the Mariner Spacecraft were instrumented at the 
spacecraft/launch vehicle interface. The instrumentation 
consisted of six (6) linear accelerometers placed on the 
adapter, such that the equivalent translational and rotational 
accelerations at the base of the spacecraft could be calculated 
assuming a rigid interface. These flight data proved extreme. y 
valuable and cost effective in the estimation of launch loads 
and test criteria for subsequent Mariner Spacecraft. 
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III. SPACE VEHICLE SYSTEM TRANSIENT LOADS ANALYSES 

Space vehicle system transient load analysis is defined as a process 
wherein a detailed model of the space vehicle system, a coupled dynamic 
model of spacecraft and launch vehicle is analytically excited by a set of 
transient forcing functions to obtain spacecraft loads. The output of 
such analyses includes detailed member loads which are used directly for 
the design, and/or acceleration of masses from which member loads can be 
derived. These analytically determined spacecraft structural loads are 
used as the basis for the structural design and qualification. This 
approach was used for the Viking Orbiter System (VOS) which was part of 
the Viking project managed by the Viking Project Office (VPO) at Langley 
Research Center (LaRC) for NASA. 

The project decision to rely almost exclusively on analysis for spacecraft 
design was a clear technological advancement in spacecraft structural 
design philosophy, and was an indication that the quality of analytical 
tools and the understanding of launch vehicle system dynamics had reached 
appreciable maturity. Further advances in computer technology, the 
development of versatile finite element programs such as NASTRAN [25], new 
methods for analyzing large dynamic systems [26], and dedicated analytic? < 
tools for booster analyses made a reliance on analysis possible. The 
approach selected for Viking was based upon the following considerations: 
(1) requirement for a lightweight structure, (2) high reliability for two 
Viking missions, (3) a new, not previously flown launch vehicle system, 
and (4) availability of launch vehicle engine forcing function data from 
previous Titan and Centaur flights. The Viking orbiter design_experience 
proved extremely valuable both from a technical and organizational point 
of view. Lessons learned from Viking were later used for developing 
simplified and more efficient methods for loads analysis. The Viking 
design process has been well documented [27], [28], [29], [30]. The 
Viking experience will be summarized here emphasizing analytical and .tost 

methods used. 
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A. Viking Orbiter Structural Design Criteria 

1. Primary Structure 

Any structural element whose failure would result directly in an 
overall spacecraft structural failure was classified as primary 
structure. Early in the program specific hardware items were 
established to be designed by transient loads analysis. 
Engineering judgement was used to select hardware with design 
loads in the low-frequency (0 to 30 Hz) range. Design loads for 
the primary structure were obtained from space vehicle system 
loads analyses not verified by test. Flight load predictions 
were obtained using mathematical models verified by a test 
program. The qualification test levels for primary structure 
were based on predicted flight loads. Wherever statistical 
information was available the design and flight loads at the 3o 
level were used. 

2. Secondary Structure 

All other structural parts whose loading was governed by the 
middle or high frequency range were designated as secondary 

structure. Secondary structure was designed to specified test 
levels. 

B. Viking Orbiter Dynamic Analyses 

As indicated in Figs. III-l and III-2, upward of nine organizations 
were responsible on the Viking project for hardware or integration 
functions which directly affected the evaluation of dynamic transient 
loads. The significant forcing functions anticipated during the 
flight are illustrated in Figures 1 1 1-3 and II 1-4. Table III-l shows 
the number of forcinq functions or conditions for which the space 
vehicle system had to be designed. A major analytical effort was 
expended on Viking to obtain design and flight loads both by the 
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Figure III 3. VOS flight loads and environments, 0-450 sec 
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Table IT.T-1 . Events for VOS Load Analysis 


No. of forcing 
functions 

Events or conditions 


Ground conditions 
Stage 0 ignition 
Airloads 

Stage 0 max acceleration 
Stage I ignition 
SRM separation 
Stage 1 burnout 
Stage II ignition 
Stage II burnout 

Centaur main engine start I (MES I) 
Centaur main engine cutoff I (MECO I) 
Centaur main engine start II 
Centaur main engine cutoff II (MECO II) 


6 

21 

5 

1 

12 

1 

29 

3 

19 

1 

1 

1 

1 
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launch vehicle contractors, Martin Marietta Aerospace (MMA) , General 
Dynamics Convair (GDC), and by JPL. The main load analysis cycles 
were performed by MMA. JPL performed simplified in-house loads 
cycles to support certain orbiter design activities using the output 
from the space vehicle system analysis. It was recognized on Viking 
from the beginning that the loads analysis process is >n iterative 
one and that as a result the design loads and hence the structural 
design would be changing between iterations. 

1. Dynamic Models 

Dynamic modeling for Viking orbiter made extensive use of modal 
synthesis and substructuring techniques 117], [20], [31]. The 
Guyan reduction technique [ 32 J was used extensively to obtain 
reduced mass and stiffness matrices. Modal truncation was 
performed based on frequency and component effective mass C33J, 
[34], Using this technique orbiter subsystems having originally 
up to 4000 static degrees of freedom each were reduced to 
approximately 300 dynamic degrees of freedom on the orbiter 
system level. A total of seven (7) subsystems constituted the 
Orbiter dynamic model C29], Throughout the dynamic modeling the 
most important parameter was the Load Transformation Matrix 
(LTM). A typical Viking orbiter dynamic model required 450 load 
transformations. Dynamic models were generated in support of 
the loads analysis cycles as well as for various subsystem and 
system test configurations. 

2. Loads Analysis Cycles 

Load analyses were performed by MMA and GDC, the launch vehicle 
contractors. Some additional analyses were performed at JPL 
using the output from the system loads analyses. 
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a. System Loads Analyses 

There were five (5) major loads cycles as shown in 
Table 1 1 1-2. The data flow for these cycles is shown in 
Figure 1II-5. Modal coupling techniques were used for 
combining the spacecraft and launch vehicle models. 

Transient responses were obtained by integration in the 
time domain. Mathematically the solution was 
straightforward, since the models were linear except for 
the Stage 0 event during which the system boundary 
conditions changed from a cantilevered model including the 
launch complex to a free-free model. Output from these 
analyses consisted of maximum/minimum data for member 
loads, accelerations and displacements .as well as time 
history data for selected parameters such as interface 
accelerations and generalized system coordinates. The 
latter data were used to perform additional analyses at 
JPL. Another exception to the linear transient response 
solution was the calculation of the structural response due 
to the transonic flight event. Member loads for this event 
were calculated by adding gust, buffet and quasi-static 
loading in a statistical manner accounting for 
dispersions l 35 1 . 

b. Simplified Spacecraft Loads Analyses 

Simplified transient spacecraft loads analyses were 
performed in-house at JPL to support Orbiter design 
iterations. Experience with such analyses showed that when 
the acceleration time history at the spacecraft/launch 
vehicle interface was assumed invariant from the old to the 
new spacecraft configuration, small changes in the 
spacecraft model sometimes produced large and unrealistic 
member loads. This was attributed to the impedance 
mismatch at the spacecraft/launch vehicle interface. To 
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Analysis with a preliminary VOS configuration 
to establish design loads 
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Analysis with a final VOS configuration to 
update design loads 

Analysis with a VOS model based on a firm 
design prior to hardware fabrication to con- 
firm design loads 

Analysis of VOS model based on modal test- 
data to establish flight loads for structural 
qualification tests 

Analysis of Viking Dynamic Simulator to ver- 
ify the load analysis process 

Analysis of the two-V-S/C flights to verify 
the load analysis process 
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Figui e 1TT-5. Dam t low lor Viking load analyses 

alleviate this problem, part of the launch vehicle was 
included in the spacecraft model and the acceleration time 
history at this interface, obtained from the previous 
analysis was used. This seemed to give much more realistic 
answers since the assumption of invariance of the 
acceleration time history had more validity as the point of 
application for the forcing function was further removed 
from the changes of the model. 

The technique of applying acceleration time histories from 
a previous analysis to a new model to obtain approximate 
member loads was called "mini-loads analysis." 

3. Error Contributions to the Load Analyses 

Since the models and the load analysis process were not perfect, 
especially for obtaining the design loads which were calculated 
from a model which had not been verified by test, a factor of 
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safety, referred to as a Loads Analysis Factor (LAF) to account 
for the uncertainties of Table 1 1 1 - 3 was used. The value for 
this factor was chosen as 1.3 based upon a qualitative 
evaluation of very limited data comparing available flight data 
with analytical predictions. The load analysis factor was 
decreased during the program as the uncertainties of Table II 1-3 
were diminished. 

C. Structural Testing 

As in all other JPL projects two types of structural testing were 
performed on the Orbiter: developmental testing and qualification 

testing. In addition the whole launch vehicle system was flight 
tested using a Viking dynamic simulator as a payload. 

1. Developmental Testing 

The objectives of developmental testing were aimed at the 
improvement of the mathematical models used for loads analyses. 

a. Modal Tests 

Extensive modal testing of all subsystems as well as the 
Viking Orbiter System was performed using the Multi-Point 
Sine (MPS) excitation technique. These tests yielded 
frequencies, mode shades, and modal damping values as a 
function of excitation levels [36l. The test results were 
correlated with analysis [37], and modifications to the 
analytical models were made. 

b. Fluid Dynamic Tests 

Since an appreciable fraction of the Viking spacecraft 
consisted of liquid propellant, the dynamic behavior of the 
fluid was a very important modelling consideration. Special 


Table 111-3. Error Contributors to VOS Load Analysis 


Spacecraft mathematical model variation 
Launch vehicle mathematical model variation 
Centaur standard shroud mathematical model variation 
Definition of the forcing functions 
Load analysis process 
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fluid dynamic tests to determine the effective mass of the 
fluid as a function of frequency of excitation and the 
ullage condition were performed [38], Analytical methods 
for obta iiing effective fluid mass data from modal tests 
were derived [39]. 

Qualification Testing 

Qualification testing consisted of a static test and a vibration 
test. 


a. Static Test 

All primary structure was qualified by a static test [40]. 

b. Vibration Test 

Secondary structure was qualified during a three axis 
sinusoidal test [41] and an acoustic test. The input to 
the sinusoidal test was limited so as not to exceed design 
loads in the primary structure. A thorough analysis of the 
test configuration was performed and excellent correlation 
of the analytical predictions with test responses was 
obtained [42] . 

Flight Testing 

A Viking Dynamic Simulator (VDS) was used as u payload to test 
the newly integrated Titan/Centaur launch vehicle system [43]. 
The objectives were to determine the accuracy of the loads 
analysis process and to verify the orbiter flight data 
instrumentation plan [44]. Excellent transient flight load data 
were obtained from the VDS program which was used to verify the 
analytical predictions. 
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D. Flight Instrumentation 

Fliqht instrumentation played a very important role in the Viking 
program. As discussed earlier the VDS flight measurements were used 
to verify the loads analysis process. Flight instrumentation on the 
Viking spacecraft consisted of strain gauges and accelerometers 
£ 45 ], ’these data were used to assess the margins of safety of the 
first Viking spacecraft before the second spacecraft was launched. 

Such response measurements were also used to synthesize launch 
vehicle forcing functions [46], The flight measurements obtained 
from VDS showed an appreciably higher loading than predicted by 
analysis for the Stage 0 Ignition event. The VDS flight measurements 
were used to re-evaluate the Stage 0 Ignition forcing function. 

E. Summary of Experiences Gained From Viking 

As discussed above the Viking spacecraft structural design was an 
integrated effort both organizationally and analytically. As a 
result, the Viking spacecraft structure was designed to loads which 
were very close to those measured in flight [47], [48]. Appendix A 
summarizes the comparison of fliqht measured data with predictions 
for various events. This integrated loads analysis approach is 
technically sound and produces a lightweight highly reliable 
structure. To achieve this several prerequisites have to be met. 

In the technical area, reliable, verified launch vehicle models and 
forcinq functions have to be available. This implies a previously 
flown launch vehicle system. On the management level this approach 
requires a commitment both in manpower and resources and the 
recognition that this is an iterative process. 

The numerous organizational and technical interfaces inherent in the 
space vehicle system transient loads analysis approach resulted in 
appreciable elapsed time between the generation of a spacecraft model 
and the availability of spacecraft leads. The output of the various 
load cycles did not always support the spacecraft design effort in a 
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timely fashion. Furthermore, it was observed that some member loads 
were sensitive to small design changes. 

These experiences prompted the search for a new methodology to obtain 
spacecraft member loads. The objectives were to simplify the loads 
effort by decoupling the procedure as much as possible from the 
launch vehicle system and to search for methods to determine an upper 
bound for the loads. This seemed especially appropriate with the 
advent of the Space Transportation System (STS), a new, and as yet 
not flown reusable launch vehicle system for which the dynamic model 
has not been confirmed and the forcing functions have not been 
measured in flight. 

The next few sections will summarize the methods being proposed. 

Some of them are currently being used to obtain design loads for 
payloads to be launched by the STS. 
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IV. LOADS PREDICTION METHODS FOR SHUTTLE PAYLOADS 
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design loads for the early STS interplanetary spacecraft. 

1. Preliminary Design Loads - The curve and mass 

acceleration curve 

For preliminary design purpose an upper bound of payload dynamic 
response as a function of payload or component effective weight 
has been estimated. These semi-empirical estimates are based 
upon the responses of various spacecraft obtained from 
experimental, analytical and flight data for a variety of 
payloads flown on various expendable launch vehicles of the past 
and from analysis done on the STS. 

The purpose of developing simplified rules for estimating 
preliminary loads is to give the analyst and the designer a 
"rule of thumb" for sizing of a preliminary structure and to 
assist load estimation in design studies. These methods are 
also intended to be used for the design of secondary structure 
for which the load analysis approach is not applicable. 

The preliminary loads are estimated on the basis of a two 
degrees of freedom system subjected to an impulse function as 
shown in Fig. IV-1, where one of the degrees of freedom 
represents the payload or component weight and the other 
represents the launch vehicle. It can be shown that the maximum 
response can be expressed as 

(IV-1) 



where C a constant, w = weight (lb) 
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Figure IV-1. Two Degrees of Freedom System Subject 
to Impulse Function 


The constant C depends on the weight of the launch vehicle. 
Extensive data from previous programs have estimated these 
constants, C = 185 for Atlas-type vehicles and C = 220 for Titan 
vehicles. Using the shuttle load factors, it was originally 
estimated that C = 285 for the shuttle payloads. More recent 
analysis shows C = 230 for the shuttle payloads. 

A disadvantage of this relationship is that as the component 
weiqht tends to zero, the response or load factor tends to 
infinity. Using the derivation of Ref. 51, it can be shown that 
m the limit, w 0, the response can be expressed as 

Limit x = i Xo i 

2.718 | 2|max 


w -* 0 


(IV-2) 
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where Q = gain factor = ^ 

a = percent critical damping, C/C c 

x 2 = interface response of launch vehicle 

carrying no payload (unloaded launch vehicle 
response) 

Using the outlined approach, a load factor curve can then be 
constructed such as the one shown in the Fig. IV-2 in which the 
preliminary upper bound response for the Galileo spacecraft 
components are plotted as functions of the effective weight. 

The upper weight range is impedance or weight controlled with 
very little effect of damping and in the lower weight range it 
is damping controlled. 

Currently, this approach is being used at JPL for preliminary 
design purposes and for estimating load factors for the design 
of secondary structures for Galileo and International Solar 
Polar Mission (ISPM) and other STS launched payloads. Figure 
IV-2 shows typical results obtained by this approach for the 
Galileo spacecraft. 

2. Generalized Shock Spectra Method 

A more accurate method for obtaining upper bound loads is the 
shock spectra method. Estimates of loads are obtained by 
combining the maximum responses for various combinations of 
launch vehicle and spacecraft normal modes, and by allowing for 
impedance effect between launch vehicle and spacecraft and 
frequency shifts to obtain the most adverse combination of 
dynamic response. The goal is to develop a cost effective tool 
for obtaining design loads in a timely fashion. The method is 
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low in cost as compared to a transient load analysis, it is 
insensitive to design changes, and hopefully moderately 
conservative in load prediction. 

The shock spectra method was originally presented in 
Reference [49], [50] and later expanded as a generalized shock 
spectra in Reference [51]. 

The following general observations are fundamental to 
understanding the rationale of the shock spectra approaches: 

(1) The general objective is to avoid the cost of a launch 
vehicle/new spacecraft overall transient analysis. 

(2) A previous dynamic analysis has been performed for the 
same launch vehicle with another or a dummy 
spacecraft, the results of which are used as inputs to 
the generalized shock spectra method. The method is 
derived to utilize modes of the launch vehicle loaded 
by a rigid mass only, or unloaded. If such an 
analysis is not available, the proper data can be 
recovered from the analysis of the launch vehicle 
loaded with a spacecraft having rigid or elastic mass 
representation. 

(3) The structural analyst needs only to determine the 
worst case load maxima or the bounds, rather than time 
histories, of the structural response. Since maxima 
or bounds are the objective of the shock spectra 
concept, a shock spectra approach is readily 
applicable to the structural design process. 

(4) Any loads analysis - transient, shock spectra, or 
other - incorporates two basic items: a model 
idealizing the dynamic environment, i.e., the forcing 
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function, and another idealizing the composite structural 
system. The form of each model is influenced by the 

selected approach, while its complexity is constrained by 
cost and time. 


(5) The original launch vehicle forcing functions, usually 
unknown to the spacecraft analyst, may be idealized by 
a simpler form to obtain an explicit closed form 
solution. A complete definition of this idealization 
is possible from the previous analysis of item (1) 
above, and is done at t.he modal level. 

Only the highlights of the generalized shock spectra method will 
be summarized here. The reader is referred to Reference 51 for 
a detailed derivation of the method. 


First, the modal forcing function F*(t) corresponding to the 
modes of the launch vehicle loaded by a rigid mass (or unloaded) 
at the spacecraft interface and representing the modal 
contribution of an actual flight event, is modeled, regardless 
of its physical point_of application, by an -equivalent launch 
vehicle modal forcing function, F (t). Unlike the actual 
complex transient force, the equivalent forcinq function assumes 
a simple form of variation with time. Here, an impulse delta 
function F g Jt) - Fg^t) with a yet-unknown magnitude 
F(D v or equivalently a velocity with magnitude v n , is 
chosen for convenience. ^ 


The choice of a simplified forcing function as an impulse 
emphasizes the view that the shape of the response time-history 
is of little consequence, and that only the peak or bound of the 
response is of interest. Therefore, any forcing function that 
would reproduce a response with the same maximum peak or bound 
as the actual forcing function is acceptable. The equivalency 
between the actual forcing function and the idealized one is 


established, not on the basis of producing identical response 
time histories, but on the basis of producing an identical peak 
of the shock spectra of each of the launch vehicle modal 
response qj(t), Figure 1 V- 3 , derived from the previously 
performed launch vehicle/dummy spacecraft analysis. An 
alternate to the direct use of q^(t) is discussed later. 

The use of modal shock spectra, rather than the interface 
degrees-of -freedom shock spectra is significant because it 
automatically accounts for the matching of all interface 
physical deqrees of freedom, and allows one to determine the 
modal maqnitude of the impulse F e ^(t) or velocity Vg^. It 
is noteworthy that the above process of establishing the 
equivalent idealized forcing function requires knowledge of the 
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modal properties of the launch vehicle with or without a 
payload. Such information is usually available from the launch 
vehicle organization. Also note that the rigid mass at the 
launch vehicle/spacecraft interface does not have to be the 
total mass of the spacecraft to be analyzed, but can have any 
value convenient for purposes of the transient analysis 
performed earlier on the launch vehicle. This approach does 
correct for whatever mass value was used previously. 

Second, in considering the composite structural system which 
consists of a spacecraft modeled by S-normal modes and a launch 
vehicle modeled by L-normal modes, there will be (S + L) modally 
coupled equations of motion. The coupling arises because the 
S-spacecraft modes arid the L-launch vehicle modes are obtained 
from two separate modal analyses, rather than from an analysis 
where the two models are integrated. Unlike the transient 
analysis where the solution is expressed in the complete (S + L) 
space of modal coordinates and time, a bound on the complete 
solution is established by: 

(1) Idealizing the totality of (S + L) mathematical space 
of modal coordinates by an array of nested (S x L) 
mathematical subspaces, in each of which only one 
spacecraft mode is individually coupled with one 
launch vehicle mode. In this fashion, each spacecraft 
mode is couDled with L-launch vehicle modes one at a 
time. To derive a bound on the total solution in the 
original (S + L) mathematical space, an explicit 
solution in the form of spacecraft modal response time 
history is first derived for the pair of modes in a 
typical subspace. The explicit form of this solution 
is based on the idealized modal forcing function just 
discussed. Furthermore, because the spacecraft member 
loads are the objective of the analysis, and since 
these are proportional to the generalized modal 
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displacements, the response quantities used here are 
generalized modal displacements, q ? (t), from which 
an expression of the bound Q s? of q ^ ( t ) is then 
derived. 

(2) Numerical computations are made first to establish the 
bounds, Q $|i , on each of the (S x L) discrete modal 
responses. Each bound Q Sii corresponds to one of the 
(S x L) subspaces. To account for unknown design 
tolerances and variations in the structural model 
idealization, worst cases can be provided by allowing 
realistic possible tuning between each spacecraft mode 
and its nearest launch vehicle mode, and by scaling 
the entire frequency spectrum of the launch vehicle 
with respect to that of the spacecraft. 

(3) Next, a bound on the total spacecraft modal response, 
Q $ , is constructed by summation over all the 
discrete L-bounds for that spacecraft mode. The 
summation can be over absolute values, or in a 
root-sum-square sense that can also be weighted to 
account for phasing. 

(4) Finally, spacecraft member loads are obtained by 
adding the contributions of all spacecraft modes, 
either in absolute value or in a root-sum-square sense. 

In the procedure outlined above, steps (1) and (2) derive 
expressions for bounds on the discrete response, Q s(i , in each 
of the (S x L) modal subspaces, while steps (3) and (4) 
construct a bound on, Q Si , , the complete solution in the total 
(S + L) space from the discrete Q s ^ bounds. In this manner, 
much computational effort is saved over the usual transient 
analysis. 
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Because the spacecraft member loads are the object of the 
analysis here, and since these are proportional to the 
generalized modal displacements Q $ , bounds on the spacecraft 
member loads are expressed in terms of the bounds, Q , of the 
generalized modal displacements for each subspace. A typical 
generalized modal displacement bound, Q sj ^, is derived here for 
a typical subspace in which one spacecraft mode q $ is coupled 
with one launch vehicle mode q^. The governing equation is in 
the following form C51J 
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where 

= eigenvalue of the launch vehicle mode q^. 

? 

u.- s = eigenvalue of the spacecraft mode q s . 

percent of critical modal damping for the launch 
vehicle mode q^. 

= percent of critical modal damping for the 
spacecraft mode q $ . 


-40- 


< ^ 4 j >s values of the launch vehicle mode q £ at the j th 
degree of freedom. 

|Fj| = time history vector of forces applied at the 
launch vehicle degrees-of-freedom. 

and jm^j, £m.j S Jare the mass and the rigid-elastic terms 
defined in [51]. 

As mentioned before, the modal forcing function <4 > £ j> 
is modeled by an equivalent impulse delta function having a 
magnitude Fg^, or alternatively, an initial velocity with a 
magnitude Vg^. 

Then the bound, Q SJ ^, of the modal response q s of a 
spacecraft mode can be obtained as 
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The definitions for m $t ,, M and D max can be found in 
Reference 51. 

During the early stages of the design, the spacecraft and launch 
vehicle modes and frequencies are obtained from analyses that 
usually contain a large degree of uncertainty. To account for 
such uncertainties, one may introduce an artificial tuning 
between the spacecraft and launch vehicle modes. Two forms of 
artificial tuning have been identified: global and local. In 

global tuning, the entire spectrum of launch vehicle frequencies 
is incrementally scaled in either direction relative to the 
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spacecraft frequency spectrum. For each increment, a global 
response 

fCIX < iv - 5 

VS £ 

can be computed and used as a measure for determining the worst 
case for design purposes. In this scheme, tuning is achieved by 
finding the amount of relative scaling that maximizes Q. 

Clearly, limits on the allowable relative scaling must be 
selected in advance, and the search for the maximum Q - conducted 
within these limits. An alternate is to evaluate member loads 
or acceleration for each tuning. 

In the local tuning, the response is maximized for each 
spacecraft mode, one at a time. This is achieved by allowing 
the nearest launch vehicle frequency to coincide with that of 
the spacecraft frequency under consideration, provided that the 
two were originally separated by no more than a preselected 
amount. 

Since each Q SJ ^ results from coupling between a spacecraft mode 
s and only one launch vehicle mode £ , and since a complete 
representation of the launch vehicle includes more than one 
mode, say L modes, contributions due to all L modes should be 
included. The upper bound of the contribution by all the launch 
vehicle modes can be obtained in the root-sum-square sense. 
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( I V-6 a) . 


or 



(IV-6b) 


where V^w s , co^) is a weighting function that can be used to 
account for time phasing between the launch vehicle modes. 

Using Eq. (IV-6a) or (IV-6b), a bound on the relative 
displ acement | D b | is found from 

N ■ VpE 1 1 

or 

V *3 


in which 

j D bs) ■ {V «s} 

t h 

1 b s} s spacecraft modal displacement at the b DOF 
4 for the s-spacecraft mode. 

W (u> s ) = a weighting function dependent upon the 
spacecraft frequency 

A bound on the absolute acceleration a fa j for a 
degree-of-f reedom, b, on the spacecraft is found from 
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where 
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S = 1, 2, . . . , S <_ S 

L^biJ = rigid body modes of the spacecraft due to unit 
deformations imposed one at a time at the 
spacecraft-launch vehicle interface i 

O^i]{ w l v 0)l} = acce l era tions at the 
interface i for each 
spacecraft mode s. 

In Eq. (IV-8), the || — || means that a bound similar to that 
of Eq. ( I V- 7 a) and (IV-7b) is taken for the enclosed 
quantities. It is also noted that the underlined term in Eq. 
(IV-8) is a correction term that accounts for errors introduced 
if a truncated set of spacecraft modes J is used in place of the 
total S. It can be shown that when all spacecraft modes are 
considered, i.e. S = S, the correction term vanishes. 

If the modal displacement method is used to calculate the member 
forces |F a | associated with the a th force component, one may 
write 




[ c *b] pb.] w 


where 


s = 1, 2, ... S 


( IV-9) 



matrix of force coefficients whose elements are 

A. L 

the a un force component associated with 
displacements in the b ^ degree-of-freedom. 
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Again, if a truncated set of spacecraft modes S < S are used to 
compute the forces the same correction term of 

Eq. (IV-8) can be introduced to give 


[ C ab] [*bs] { Q s} 


Mps ]" 1 m-m mn 


(IV-10) 


As before, the underlined term is the member load correction due 
to mode truncations, and || — || means that a bound is taken 
during summation over the spacecraft similar to Eq. (IV-7). 
Member loads computed by the generalized modal shock spectra 
method are, by intent more conservative than those calculated by 
the transient method. 


A comparison of acceleration levels obtained by the generalized 
shock spectra method to those from the mass acceleration curve 
are contained in Appendix C. 


Transient Load Analysis Techniques 

For more accurate loads prediction, the launch vehicle/payload 
composite model will be subjected to appropriate external forcing 
functions for the response and member load calculations. These 
external forcing functions are derived from the relevant dynamic 
environments representing the launch vehicle thrusts, staging events, 
aerodynamic loads and others. Since the mathematical model is based 
on the design of the structures, this design/ analysis approach is an 
iterative process. This design/analysis process involves the inte- 
gration of the payload model and launch vehicle model and the 
subsequent transient analysis. Considerable time and cost are 
required in the integration of the composite mode! and the response 
analysis, since typically the payloads and the launch vehicle are 
developed and designed by separate organizations with their res- 
pective mathematical models which involve different computer codes, 
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coordinate systems, units and normalization procedures. In a 
previous experience, the Viking project, as many as ten organizations 
were involved in generating the composite model and its integration. 
One design/analysis cycle required up to six months [28] . 

io reduce cost and schedule, several approximation techniques have 
been developed in the past. One method utilizes the fact that 
payload design changes during the design/analysis iteration are 
"small", such that a cost effective perturbation procedure can be 
applied to update the response analysis due to design changes [52]. 
This method was further extended to estimate the modifications of 
launch vehicle/payload interface responses due to payload design 
changes [53]. The modified interface responses can then be applied 
to the base of the payload to calculate the payload responses and 
loads. This method is not only cost effective but also can be 
implemented by the payload organization without interfacing with the 
launch vehicle organization. However, this method is developed with 
the assumption that the payload is always much "smaller" than the 
launch vehicle. For shuttle-launched payloads, the effects of 
payload dynamics can no longer be considered as "small" for future 
dynamic events such as the abort landing of a fully loaded shuttle 
orbiter. in these events, the interaction between the shuttle 
orbiter and its payload is critical to both the orbiter and payload 
loads [54] . Therefore it seems that the launch vehicle/payload 

composite model integration and its transient analysis cannot be 
avoided. 

The simplified loads analysis methods aim at simulating a time 
response solution within the payload organization with minimal launch 
vehicle information. The objective is to greatly reduce the 
interface between the payload and launch vehicle organizatuon. Thus 
methods to obtain the payload response of a launch vehicle/payload 
composite system from the results of another launch vehicle/payload 
composite system under the identical forcing function have been 
developed. Fig. IV-4 shows two composite systems with identical 
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where 


I* 1 !’ jy^ = launch vehicle degrees-of-freedom (DOF) of the 
compxisj±e-system (I) and (II), respectively. 
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M. (y 2 ) ■ payload OOF of the composite system („ and 
UU, respectively. 

[■")] ■ mass matrix of the launch vehicle. 

[*=,] ■ stiffness matrix of the launch vehicle. 

fc]. ft] ■ «ss matrix of the payload a and payload B of 
e system (I) and (II), respectively. 


[ k n]’ [ k l2]> [ k 21 ]. [ k 22 ], 

N- N- [*„]. [y. 


F(t) = 


sub-matrices of the total payloa< 
stiffness partitioned into launcl 
vehicle/ payload interface DOF an< 
payload DOF for the system (I) ar 
(II), respectively. 

vector representing the external 
farcing functions acting on the 
launch vehicle DOF. 


Assuming a “partial" solution of Eg. (IV-11) is available, the 

^•"v-l^ytinrth: " 

>e implemented by the payload organising alone " Pr ° Ce ” " 

Although damping is not included in Egs. (IV-11) and (IV-12, , t 

slplicUy^iT 'n'p ‘ he f °™ ° f "° da1 dampin9 ' fl,s ° f0l “ 
y. will be assumed at present that the payload is 

supported ,n a statically^^,,^ manner such that 
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[ k n] ■ h ] 1 |>] M 

[ k 2l] ’ • [♦r] 7 [ k 22] " [ k 12] T 

where 


(IV-13) 


R, 1 = payload A rigid body transformation matrix defined as 
L Rj the payload displacements due to unit displacement of 

the launch vehicle/payload interface DOF, j Xjj . 

I x \ * launch vehicle/payload interface DOF connecting 
' payload to launch vehicle, a subset of the launch 

vehicle DOF | x^j . 

Next, the motion of the payload will be decomposed into two parts, 
namely, the rigid body motion and the elastic motion: 


{ x 2 } = [>R] ( X l) + ( X e ; j 


( IV- 1 4 ) 


The first term on the right-hand side of Eq. (IV-14) is the 
rigid-body motion. The second term |x e J is the elastic motion, or 
relative motion with reference to the interface. It should be noted 
that only the elastic motion (xl will generate internal loads in 
the structure. Using Eqs. ( IV-13) and (IV-14), Eq. (IV-11) can be 
transformed into the following form 


'Wr W 


X 1 

+ 

'k t 0 


x i 

Ij 

F(t) 

J"2*R m 2 


x e 


0 t<22 


x e 

| 1 

0 

\ 


where 

[v] =N t NM 

denoted as rigid-body mass of payload A. 
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Simi 1 arly, Eq. (IV-12) can be transformed into the following form: 




Vr 


where 



(VI-16) 


M ~ P a y load ® ri g id- body transformation matrix defined 
similarly as for payload A. 


[ S rr]‘C*R] T [" 2 ]Cf R ] 

denoted as rigid-body mass of payload B. 

Also the motion of the payload B is decomposed into two parts similar 
to that of payload A, Eq. (IV-14), as: 


W'tylM+kl 


(VI-17) 


where | yj j is a subset of y^j defined as the launch 
vehicle/payload interface DOF similar to that of 


Using the basic mathematical formulations presented above, two 
different approaches for obtaining transient spacecraft responses 
using simplified methods will be discussed. 


1. The Rigid Body Interface Acceleration Method 


Since a significant coupling between the shuttle orbiter and the 
payloads exists, the launch vehicle organization will usually 
perform a dynamic analysis of the shuttle/payload system with a 
rigid payload model for the purpose of verifying the launch 
vehicle integrity. A method has been developed by which the 
launch vehicle/payload interface responses are modified such 
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that the effect of payload elasticities .are taken into 
consideration [55]. The payload responses and loads can then be 
obtained by applying the modified interface responses to the 
base of the payload. The advantage is that the entire procedure 
can be implemented within the payload organization such that 
timely design/analysis iteration can be performed. However, the 
disadvantage is that only the analytical interface responses can 
be used, since it is unlikely that a rigid payload will be 
flown. Therefore, the measured flight data cannot be directly 
applied in the design/analysis process. 

For a rigid payload, it is postulated that the stiffness of the 
payload is infinitely large, such that no elastic motion of the 
payload can be realized, i.e., | x e | = 0. Then the governing 
equation (IVJL5) can be written as 


_”l + m rr] { X l} + [ k l] M ' kk < VI - 18) 

Physically, [m ] represents the distribution of rigid payload 
mass onto the interface degrees-of-freedom. Generally, for a 
typical payload, non-structural weights such as instrumentation, 
electronics, propellants, etc., constitute the major portion of 
the payload mass and the weight of the load carrying structure 
is only a small portion of the total payload mass. Therefore, 
early in the project the payload rigid-body mass [m rr ] can be 
estimated prior to the actual design since only the mass 
distribution and geometric configuration are required to 
establish pn rr ] and [<t> R ]. It is a common practice that the 
payload organization will provide an estimated [m rr ] to the 
launch vehicle organization early in the project such that the 
launch vehicle/rigid body payload composite model, i.e., 

Eq. ( I V- 18 ) , can be constructed. The main purpose of such a 
model is the verification of the launch vehicle loading. 
Meanwhile the interface accelerations of the model can be 
obtained. 
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The governing equation for the launch vehicle/rigid payload 
composite model, Eq. (IV-18) can be reduced to the generalized 
coordinate formulation as: 


+ 2 


S J['“u] {*il + ; 


where 




(VI-19) 


IM ■ W t v i) ) 

> (IV-20) 

(G(t)j - (t) J ) 



m rr] [*l] = “ 



unity matrix 


(IV-21) 


Clearly, and['w^ ] are the eigenvectors and eigenvalues of 
the launch vehicle/rigid payload composite model, respectively. 
Also it should be noted that modal damping has been included in 
Eq. (IV-19) and the elements in ^ J represent the 
percentage of critical damping for each mode. 

It is important that the launch vehicle organization has 
obtained the solutions to the launch vehicle/rigid payload 
composite model, Eq. (IV-18), and made them available to the 
payload organization. Precisely, the following are required: 

:he interface accelerations of the events under consideration, 

Xj J, the eigenvalues of the launch vehicle/rigid payload 
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composite system , the corresponding eigenvectors of the 
interface degrees-of-freedom, jjJ'j [^ 2 ! a subset of 

and the modal damping £p 1v . The objective of the 
method is to use these data provided by the launch vehicle 
organization to calculate the elastic payload responses and 
loads under the same dynamic events, without having to solve a 
new launch vehicle/elastic payload composite model similar to 
Eq. ( IV-16) . 

First a transformation will be defined as 


i*n pi 0 

U) ’[° *2 

where £<t>^ are the eigenvectors of the launch vehicle/rigid 
payload system, Eq. (IV-21), and [d> 2 ] are the eigenvectors of 
the elastic payload constrained at interface. ^ 2 ^ satis ^ 
the following orthogonality conditions: 


n 

(u 2 ) 


(IV-22) 





unity matrix 



eigenvalues of elastic payload 


( I V- 23 ) 


The governing Eq. (IV-16) can then be written as 



(IV-24J 
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where 


[■„] ■ [4 ■ Kf [V 



s p 2 55 pa.' 


payload 


modal damping 


With the information provided by the launch vehicle 
organization, Eq. (IV-24) can be constructed by the payload 
organization except the generalized forcing G(t). To eliminate 
this difficulty, a modal response due to the launch 
vehicle/rigid payload interface acceleration will be defined as 



m 


er 


■ -M H i’ll 

(IV-25) 


It should be noted that {v 2 | can be obtained once the 
interface accelerations jxjj. are available. Eqs. (IV-19) and 
(IV-25) can now be combined as 
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( IV-26) 

Then the solution to Eq. (IV-24) will be decomposed into two 
parts as 
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Upon substitution, the following equation can be obtained 
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(IV-28) 


In Eq. (IV-28), the generalized forcing function G(t) has been 
replaced by V 2 (t) which in turn can be obtained from 
Eq. (IV-25) once the launch vehicle/rigid payload interface 
accelerations are available. Thus the payload responses 
y g can be obtained by the payload organization once 
Eqs. (IV-25) and (IV-28), .are solved as follows 


W = hlKI* p] (W + W) 


The above describes a method by which a transient loads analysis 
can be performed within the payload organization. The method 
requires certain information from the launch vehicle 
organization including the launch vehicle/rigid payload 
composite frequencies, modal damping modal displacements for the 
interface degrees-of-freedom, and interface responses due to 
various events. The advantages of the method are the ability of 
the payload organization to perform complete a design/load 
analysis cycle independently, thus eliminating the costly and 
time consuming interfaces between the launch vehicle and payload 
organizations. The method developed above has no limitation 
regarding the relative size of the payload to launch vehicle or 
the type of forcing function used in the transient analysis. 

The detailed mathematical derivation and the application of the 
method to a realistic problem can be found in Reference 55. 


2. Recovered Transient Analysis 

The objective of this method is to develop a procedure by which 
the interface accelerations of one launch vehicle/payload 
system, either analytically obtained or in measured flight, can 
oe used directly in another launch vehicle/payload system design 
procedure, thus, the name "Recovered Transient Analysis." This 
effort should be especially beneficial for future payload 
designs using the shuttle as a launch vehicle. The payload 
dynamic loads can be obtained by performing a transient analysis 
of the payload model using the interface accelerations modified 
from another composite system, which consists of an identical 
launch vehicle, the shuttle, and a different payload, as forcing 
functions. In principle the proposed transient analysis 
recovers the interface accelerations for the unloaded launch 
vehicle and then modifies them to include the dynamic 
characteristics of the new payload. After a series of shuttle 
launches, the flight measured interface accelerations can be 
used to establish a payload forcing function data base for 
subsequent payload designs. The payload organization using the 
proposed loads analysis method can then perform payload loads 
analysis within the payload organization in a cost effective and 
timely manner. 

In a recent article, C563 , it was pointed out that the current 
dynamic analysis methods are unresponsive to the design process 
because of the long computational times and their associated 
costs. Their use is virtually precluded in the early design 
stages where frequent design changes are occurring. This method 
will not only simplify the analytical process thus reducing the 
cost, but also will provide results in a timely manner due to 
minimal interaction between the various organizations. 

Therefore, the rigor and potential accuracy of a systematic 
analytical procedure can be brought to bear on not only the 
design process in the early stages but also as a guide for the 
qualification testing in the later stages of the project. 
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Since the objective is to utilize the results of the Composite 
System I to solve for the response of the Composite System II, 
as shown in the Fig. IV-4, the solution of Eq. (IV-15) which is 
the governing equation for the System I will be examined. 




where [<)>] and [<t> e ] are the eigenvectors of the Eq. (IV-15) 
for the launch vehicle degrees of freedom (DOF) and payload 
elastic DOF, respectively, and |q(t)J is the generalized 
coordinate vector. The eigenvectors satisfy the following 
orthogonality conditions: 


* I I m i + m rr ^R m 2 I I * 


^e m 2^F 


> 1 T fk, 


4> e 0 


0 1 [> 


k 22 ^e 


where is_the composite System I eigenvalue matrix. 

From Eq. (IV-14), the complete eigenvectors for Payload A can be 


obtained as 


M ■ M h] + [y 


( IV-32) 
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Where M are the eigenvectors for the Interface DOF, a 
subset of the eigenvectors for the launch vehicle DOF [*J. 


By substituting Eq. (IV-30) into Eq. (IV-15) and using the 
orthogonality conditions as expressed in Eq. (IV-31), one 

obtains the governing equations for the generalized coordinates 
q(t) . 


( |q) (q j = C 4. ] T { F(t)J (IV-33) 


Note that modal damping, C'2pu)J, has been included. 

One other quantity, [m er ], representing the inertia coupling 
effects between the launch vehicle and the payload is available 
from the analysis. It is defined as follows: 


[ 5 er] ■ [fj ["2] [*r] 


(IV-34) 


summary, the following results from Composite Systei 
have to be made available to the payload organization. 


1 ) 

2 ) 

3) 

4) 


The system eigenvalues^' w 2 v j . 

The system eigenvectors at interface DOF [ 4 ^]. 

Rigid-body mass and inertia coupling matrices, 
and for]- 



ime history of the generalized coordinates |q(t)| or 
q(t)j. 1 1 
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With this information, the responses of System II represented by 
Eq. (IV-16) will be sought. First, the external forcing 
functions |F(t)| must be eliminated by combining Eqs. (IV-15) 
and (IV-16). 
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(IV-35 ) 

where 


’ Fi - x i 


)• [*r] (*i| ’ [*r] (*i 


(IV-36) 


F(t) | has been 


In Eq. (IV-35), the external forcing function ^ x 
replaced by the modal acceleration q(t)J which will be 
available to the payload organization. However, to solve for 
Eq. (IV-35), the eigenvalues and eigenvectors of System I are 
required. From Eq. (IV-35), it appears that the launch vehicle 
mass and stiffness matrices, [m-jJ and [k-^, are needed which 
means the launch vehicle model is needed. Since it is an a 
priori assumption that the launch vehicle model will not be 
available, a method will be devised by which the required 
eigendata can be extracted without the use of the detailed 
launch vehicle model. Only the generalized, or modal launch 
vehicle model will be required. Basically, the method will 
first obtain the eigendata for an unloaded launch vehicle by 
removing Payload A from the Composite System I. Next, the 
unloaded launch vehicle will be coupled with Payload B to obtain 
the eigendata for Composite System II. For reasons of matrix 
size compatibility, a system consisting of an unloaded launch 
vehicle and a cantilevered Payload A will be studied. 



(IV-40) 
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All the matrices in Eq. (IV-40) are available to the payload 
organization and they are in much simpler form than the launch 
vehicle mass and stiffness matrices. The eigenvalues and 
eigenvectors for Eq. (IV-40) will be denoted as and 

[jlj* respectively. From Eq. ( I V- 38 ) the eigenvectors for the 
unloaded launch vehicle interface DOF will be 



( I V -41 ) 


As mentioned before, the eigenvalues and the eigenvectors 
obtained from Eq. (IV-40) contain the effect of the cantilevered 
Payload A. It must be identified and removed from the results. 
One way is to use the orthogonal ity condition to identify the 
payload modes. 

Let 



( IV-42) 


From Egs. (IV-38) and (IV-42), it is clear that the cantilevered 
payload modes are among the modes Since the payload 

modes are orthogonal modes with respect to the mass matrix, the 
following multiplication will ideally produce unit diagonal 
terms and zero off diagonal terms for those payload modes. 


H ■ [vf H 


K os 


(IV-43) 


Although truncation errors will contaminate the unit and zero 
terms, the errors should be small if sufficient modes are 
accounted for. Thus the matrix can he used to sort out 
the payload modes. 


After the eigendata a" unloaded launch vehicle are obtained, 
the eigendata for the new composite model, Eq, (IV-35) 
consisting of the identical launch vehicle and a .new payload B 
can be obtained. 


'• L l 


0 1 u. 


*2 u 2 


(IV-44 


where L<J>qD and [> 2 ] are the eigenvectors of the unloaded 
launch vehicle and the cantilevered Payload B, respectively. 
Note that the matrix {^qi] is available from Eq. (IV-41). 
Therefore the following orthogonality conditions are satisfied: 


W m w ■ r ■ -j 

[*.]' M W NJ 
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[*2] _ m 2_ = 

[*2] , k 22_ 3. r*- “2 ."j 


( IV-46 


where J and J are the eigenvalues of the unloaded 
launch vehicle and the cantilevered Payload 3, respectively. 
Using the transformation of Eq. (IV-44) the eigenproblem for the 
new composite model, Eq. (IV-35) can be rewritten as 
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where 
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_m 2 ] [y 


( IV-48) 


The elements in the coefficient matrices of Eq. (IV-47) are 
typically available either from the unloaded launch vehicle 
eigendata such as the [<|> 0I ] and J and from the 
cantilevered new payload eigendata such as the [m 1 r| 1 
and L wg J. Therefore, the eigenvectors and eigenvalues can 
be calculated by the payload organization. They satisfy the 
following orthogonality conditions. 



where [^], [<|/ p ] are the eigenvectors and J are the 
eigenvalues. The eigenvectors of the launch vehicle/payload 
interface DOF are given as follows: 


.* 1 ] = *oi. h_ 


( IV-50) 


The responses of the composite model, Eq. (IV-35) will be solved 
by using the following transformation: 



0 1 1 > 


*2 K 


( IV-51 


Then the uncoupled modal equations can be obtained as 


+ 2p( 


q | + j q j - [G] jq(t)J ( IV-52 


where 


G • *, V 
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( IV-53 


Since the quantities on the right hand side of Eq. (IV-52) are 
available to the payload organization, the modal response q(t) 
can be obtained from the LV organization in a timely manner. It 
should be noted that modal damping ^2pdJ^ has been included in 
Eq. (IV-52). 

The discrete payload accelecations of each DOF can be obtained 
by using Eqs. (IV-17), (IV-36) and (IV-51) as follows: 


>] (M + IM 

[h] ( ( 2 I } + l 5 l}) * {M 


|>J M i q(t) i * 1*1 U * N |>J f q(t) | 

( kl k,l [*.l + kl fv„l 1 !q(t) ! + is,j -0V-5-I 
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The structural member loads, |pf, can be obtained from the 
following operations: 


{pf s [s] jy e j = [s] 





( IV-55) 


where [$] is the loads coefficient matrix which relates the 
elastic deformation to the member loads. Typically, in a 
payload dynamic analysis, the product [Sl^J, rather that [S3 
is obtained. 

The described recovered transient analysis should theoretically 
provide the same results as that of the full scale transient 
analysis if the damping is absent. In view of the fact that 
generally only very light dampings are involved with the payload 
structural systems, good accuracy for the results can be 
expected. Another source of error is the modal truncation. One 
should expect similar errors due to modal truncation just as any 
other dynamical system solved by the modal method. This method 
has been demonstrated by a realistic sample problem which can be 
found in Appendix D. 


C. Comparison of Load Estimation Methods to Transient Loads Methods 

Spacecraft members loads obtained by estimation methods, such as the 
generalized shock spectra/impedance method are by intent more 
conservative than those obtained by transient load analysis. It is, 
however, difficult to quantify this conservatism analytically. One 
of the reasons for this difficulty is that in order to fully assess 
the conservatism a particular spacecraft structure would have to be 
designed by both methods; this implies several design and load 
analysis iterations. Such data is not available. Some limited data 
for realistic payload structural systems from past programs is 
available and will be discussed here to assess the— degree of 
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conservatism of the load estimation methods. These comparisons 
should serve as an aid in the selection of loads methods for future 
projects. 

The earlier version of the shock spectra method [49] was applied to 
the Viking Orbiter primary structure. The loads were then compared 
with the Viking design loads which were calculated by the full scale 
transient analysis on the system model. The results have been 
reported previously [57], In evaluating this comparison it should 
again be remembered that the data presented is based on one load 
cycle using the shock spectra method for a spacecraft model designed 
by a system transient loads analysis. Thus this is not a 
comprehensive comparison of a structural system designed by the two 
methods. 

Another type of comparative data is available from Voyager. This 
spacecraft structure was designed using the earlier version of the 
shock spectra method. The design loads were verified by transient 
analyses and found to be conservative. Flight instrumentation on 
b.oth Voyager spacecraft was used to reconstruct member flight loads. 

A detailed comparison of the design loads, verification loads and 
flight loads is contained in Reference 58. 

A summary of the results of these two comparisons are shown in Tables 
IV-1 and VI-2. Structural members used for this comparison are 
primary truss elements, subjected to significant loading during the 
mission. Table IV-1 contains ratios of shock spectra loads to 
transient member loads for both Viking and Voyager primary truss 
members for two launch vehicle events, namely launch (or Stage 0 
Ignition) and Stage I Burnout. This comparison shows that loads 
calculated by the early shock spectra method when compared for each 
event can be considerably more conservative than those obtained from 
transient analysis. For individual members the conservatism ranges 
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from a favorable 26 percent to a very high 500 percent. The latter 
number is due to the assumed deleterious tuning between the launch 
vehicle and the spacecraft normal modes. A more meaningful 
comparison can be made by considering the design loads for the 
various members. Ratios of shock spectra design loads to transient 
design loads for both Viking and Voyager are shown in Table IV-2. 

The design loads are more meaningful since thay are more directly 
related to structural weight. The member design loads are obtained 
by considering all flight events, some of which do not lend 
themselves to a transient loads analysis such as the maximum aq event 
which is treated as a steady-state excitation. 

For the Viking spacecraft, an estimate of the increase in the 
structural weight was made [57] using the design loads obtained by 
the shock spectra method. A minimum additional structural weight of 
6.9 percent was estimated by using the actual Viking orbiter primary 
structure capability determined by the test program, not the 
analysis, since this was the only available data. 

In some cases, the tested capability was much greater than the 
analytically predicted values. If analytically predicted 
capabilities were used as a reference, as would be done in designing 
the structure, the weight increase would be larger. Since the 
conservatism in structural design is dictated by factors such as 
manufacturing considerations, handling, and design load conditions 
that changed in the course of the project, a weight increase of up to 
50 percent might be possible. These considerations, as they arose on 
Viking, are reflected in an overall structural weight increase of 
only 6.9 percent using a single loads cycle and the shock spectra/ 
impedance method rather than the transient loads analysis method. 

Thus the values shown in Tables IV-1 and IV-2 do not automatically 
result in a very heavy structure because: 

1. The design process does not use all the computed loads but 
only a subset of the highest loads. 


Member 

Launch 

Stage I BO 

750 

2.66 

2.70 

751 

3.64 

2.15 

752 

4.65 

1.66 

753 

3.53 

1.90 

754 

2.05 

2.25 

755 

2.98 

2.38 


Member 

Launch 

Stage 1 BO 

6801 

1.82 

1.55 

6802 

1.88 

1.43 

6803 

2.84 

2.93 

6804 

3.21 

3.62 

6805 

1.99 

1.47 

6806 

1.88 

1.26 

6807 

3.40 

5.66 

6808 

2.85 

3.27 


Table IV-2, Comparison of Shock Spectra and Transient 
Design Loads for Major Structural Elements 


Viking VLCA Truss 


Voyager Mission Module Truss 


Ratio 

Member Shock Spectra/Design Load 


Ratio 

Member Shock Spectra/Design Load 


750 

2.13 

6801 

1.82 

751 

2.23 

6802 

1.88 

752 

3.38 

6803 

1.82 

753 

3.14 

6804 

2.46 

754 

1.50 

6805 

1.99 

755 

2.52 

6806 

1.88 



6807 

2.37 


6808 


2.48 


\ 
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2. T hc shock spectra loads tend to equalize the loads in like 
members, make the loads less sensitive to directional 
effect and allow for deleterious tuning between launch 
vehicle modes and spacecraft modes. 


3 . 


The advantage of the transient analysis can rarely be fully 
exploited in design since tailoring member size to loads 
and degree of freedom, member by member, is not feasible in 
practice because of schedule and cost. In addition, design 
loads derived from transient analysis must be chosen higher 
than transient loads prediction in order to account for the 
deleterious tuning mentioned above. 


herefore, although the conservatism can be high the shock spectra 
method was found to be a cost effective loads approach. The data 
does not correspond to the shock spectra derivation of Section IV-A. 
The latter is currently being used for the design of the Galileo and 
the International Solar Polar Mission (ISPM) spacecraft. Table IV-3 
shows the comparison of the ISPM spacecraft generalized shock spectra 
loads with those obtained by the recovered transient analysis. The 
comparison indicates that the loads are indeed very closed. 
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Table TV-3. Comparison of Generalized Shock Spectra and Recovered 
Transient Loads for I.SPM 


Member 

Upper Bound 
Load 
(N) 

Maximum Load 
Prom Recovered 
Transient Analysis 
(N) 

Ratio 

Upper Bound/ 
Transient 

701 

4020 

2861 

1.41 

702 

10040 

7968 

1.26 

703 . 

12200 

10943 

1.12 

704 

7100 

7692 

0.92 

705 

4380 

3759 

..17 

706 

10260 

9430 

1.09 

707 

11500 

12197 

0.94 


703 


7100 


6891 


1.03 
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OTHER METHODS 

Other aerospace organizations are actively involved in shuttle payload 
activities. In this section four different loads methods proposed by 
various companies will be briefly reviewed. 


A. Time Uncorrelated Maxima and Minima 

The Boeing Aerospace Company (BAC) has developed a simplified method 
[59] which aims at reducing the number of normal modes to be used in 
the evaluation of the time history of the member forces obtained from 
a transient analysis. For any event, the approach is, to first 
search for the absolute maxima of each coupled generalized coordinate 
q(t). Second, the loads transformation row matrix (LTM) for each 
component is multiplied by the column of these maxima and the 
contribution to the total summation is noted and ordered by absolute 
maxima. Based on this ordering only those modes which are deemed 
important are selected. The selection criteria used can be varied 
but, in practice, those modes contributing less than a predetermined 
value, such as 0.1%, to the total summation are discarded. Based on 
these criteria typically 50% of the modes are declared trivial. The 
computed time history for the force (or moment) component is now 
calculated using only the retained modes. It should be noted that 
each force component has a different set of retained modes. This 
method saves time and seems to be an interesting innovation for the 
time domain solution technique. There is as yet a lack of data to 
fully assess the conservatism (or lack thereof) of this method. It 
should not be difficult to estimate the bound of the error. BAC 
experience indicates that for a cutoff criteria of discarding the 
0.1% contributory modes, a saving in computation time of 20 or 30% is 
achieved. For a 1% cutoff criteria the saving is approximately 75% 
in computation time. 
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B. Impedance Technique for Determining Low Frequency Payload Environments 

Martin Marietta Aerospace (MMA) has developed an "impedance" method 
[60] which eliminates the necessity of creating detailed coupled 
models as well as decreasing the scope of an overall integration 
task. This approach corrects the response of the launch 
vehicle/payload interface to reflect feedback changes associated with 
changes of the payload. All calculations are made in the frequency 
domain. The approach eliminates the necessity of computing the final 
coupled eigensolutions. The final equations are reduced to simple 
complex transfer functions. Furthermore, the launch vehicle dynamic 
characteristics required to compute these transfer functions consist 
of unloaded interface free-free modal data. By obtaining a 
"standard" set of launch vehicle models and input environmental data, 
the payload organizations should be able to calculate the expected 
low frequency environments at the launch vehicle/payload interface. 
This approach also reduces a large portion of the overall integration 
task. This approach was applied to the Long Duration Exposure 
Facility (LDEF) loads calculation for the STS environments. From the 
results, it was concluded that the frequency domain analysis was 
quite difficult to implement because of, among other things, the 
sensitivity to the damping. A comparison with the results of the 
time domain ana^sis also indicated that more meaningful information 
can be gained from the transient analysis in time domain. 

C. Coupled Base Motion Response Analysis of Payload Structural Systems 

The Marshall Space Flight Center (MSFC) has developed a simplified 
analytical approach for payload loads analysis [61]. This approach 
is similar in many respects to the recovered transient analysis 
described in Section IV B.2. The main advantage of this approach is 
the avoidance of the reconstruction of the entire launch 
vehicle/payload composite model for each loads analysis cycle. This 
conceptual approach was further studied by other investigators using 
very simple examples [62]. It was concluded that the method gave 
good results even, with modal truncation. 
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D. Modification of Flight Vehicle Vibration Modes to Account for Design 
Changes. 

As mentioned before, quite commonly, subsequent to the completion of 
a load cycle analysis, design modifications are made to the payload. 

A new load cycle is then required to assess the effects of these 
design changes on payload responses and loads. A technique has been 
developed by the Lockheed Missiles and Space Company (LMSC) that will 
bypass a portion of the new analysis cycle. The method is called 
mode modification [63], which consists of expressing the new 
composite modes in terms of the set of previous modes plus a set of 
constraint modes associated with the area of the composite system 
which is to be modified. An eigenvalue problem for the new modes is 
formulated within which the finite element mass and stiffness 
matrices in the area of modification are explicitly displayed and may 
be modified in any arbitrary manner. This method is not a 
perturbation method. Solution of the eigenvalue problem yields a set 
of new system modes of the modified composite system directly. This 
method was demonstrated using three sample problems with favorable 
results comparing them to the exact solutions. 
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VI. CONCLUDING REMARKS 

The state-of-the-art of load prediction for payloads launched by STS is 
presently relying on very crude quasi-static load factors for preliminary 
analyses and complete STS/Payload or STS/IUS/Payload system transient load 
analyses for design and verification loads. The latter process is 
expensive and time consuming due to the complexity of the dynamic models 
and the many organizational interfaces. At the present time the accuracy 
of these analyses is questionable due to the lack of definition of the 
forcing functions and uncertainties in the model of a launch vehicle which 
has not yet been flown. 

Over the last few years the problem of reducing the time and expenditure 
for the design/analysis cycle of payloads has been addressed. Design 
loads for the Voyager spacecraft were obtained in-house at JPL using upper 
bound loads obtained by the Shock Spectra/Impedance Method. The design 
loads were verified by a system transient analysis later in the program. 

A comparison of Voyager flight data with design data has shown that the 
design was more conservative than Viking, which was designed using 
transient analysis, at a substantial cosJLxeduction and some increase in 
structural weight. 

Recently the shock spectra method has been improved and it is currently 
being used at JPL to design the Galileo (GLL) and the International Solar 
Polar Mission (IS PM ) spacecraft. Designing spacecraft to be flown on STS 
to upper bound levels is appropriate due to the uncertainties of a yet not 
flown launch vehicle. 

Recognizing the necessity of simplifying the transient loads analysis 
process JPL is also developing methods for decoupling the payload loads 
analysis process from the launch vehicle as much as possible. These 
methods are aimed at reducing the cost and improving the schedule for 
transient load analyses. Ultimately, as flight data for the STS will 
become available it is expected that these simplified transient loads 
analysis methods will provide the most accurate loads prediction at 
reasonable cost. 
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APPENDIX A 


Comparison of Measured and Analytically 
Predicted Member Forces 
for Viking Orbiter 


oth viking spacecraft. Viking 1 and Viking 2, had strain gauges Installed on 
he Viking Lander Capsule Adapter (VLCA) truss. These strain gauges measure 
.he axial forces in each of the six -Mrs of the truss. 1 * *-1 «« 

■caparison of the flight data to pre-flight predictions. The fma flight 
load prediction for Viking was made using the spacecraft structural 
verified by a test (Model VIII). Only those launch vehicle events onside 
critical were analyzed using this model. Other events, considered 
critical were analyzed only early in the program using spacecra 

models (Model I and IV). The predictions for the latter Models 
expected to be as good as for Model VIII. 
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Table A- 


Membc; r/Mca s . No . 


750/CY 186S 
751/CY 187S 
752/CY 188S 
753/CY 189S 
754/CY 1.90S 
755/CY 191S 


750/CY 186S 
751/CY 187S 
752/CY 188S 
753/CY 189S 
754/CY 190S 
755/CY~191S 


750/CY 186S 
751/CY 187S 
752/CY 188S 
753/CY 189S 
754/CY 190S 
755/CY 191S 


750/CY 186S 
751/CY 187S 
752/CY 188S 
753/CY 189S 
754/CY 190S 
755/CY 191S 


750/CY 1 86S 
751/CY 187S 
752/CY 188S 
753/CY 189S 
754/CY 190S 
755/CY 1 9 IS 


1. Comparison of Vi kins 1 and Viking 2 Measured VI, CA 
Forces to Prof. light Analytica-1 Predictions 


(Compr 

Minimum 
ession) f 

lb. 


Max i mum 
(Tension), lb. 


Predicted 

Viking 1 

Viking 2 

Predicted Vildna 1 

Viking 2 

S Uige 

0 ignition (Model VIII) 


-2900 

-2000 

-2300 

900 

800 

100 

-2700 

-1600 

-1300 

1800 

800 

1000 

-2300 

-2000 

-2200 

400 

100 

0 

-2900 

-2500 

-2800 

900 

200 

1100 

-2800 

-2200 

-1500 

1900 

600 

100 

-2800 

-2200 

-zaoo 

800 

200 

200 


(Max « q 

(Model VIII) 




-3200 

-2000 

-1900 

1000 

- 500 

- 400 

-2900 

-1500 

-1200 

1800 

300 

200 

-3400 

-2100 

-2100 

1200 

- 400 

- 300 

-3600 

-1900 

-2000 

1400 

- 300 

- 300 

-3000 

-1500 

-1400 

2000 

200 

400 

-3400 

-1900 

-2200 

1200 

- 300 

- 400 

Stage I Ignition (Model 

J.) 




-1600 

-1400 

-1800 

0 

- 800 

- 800 

-1200 

-1200 

-1200 

0 

0 

0 

-2000 

-2000 

-1700 

0 

- 600 

- 600 

-2100 

-1800 

-1800 

0 

- 600 

- 600 

-1800 

-1200 

-1200 

0 

0 

- 100 

-2000 

-1900 

-1900 

0 

- 500 

- 700 


SRM -Jettison 

iJ^lodel I) 




-1200 

-1500 

-1600 

0 

- 400 

- 300 

- 700 

-1000 

-1000 

0 

200 

300 

-1500 

-1700 

-1700 

0 

- 500 

- 200 

-1200 

-1400 

-1600 

0 

- 100 

0 

- 800 

- 900 

-1.100 

0 

100 

200 

-1400 

-1500 

-1500 

0 

- 100 

- 300 


Stage 1 Burn 

(Model VIII) 





-4900 

-2200 

-2200 

1000 

-1600 

-1600 

-2500 

-1100 

-1000 

500 

- 800 

- 800 

-5400 

-2200 

-2200 

1400 

-1800 

-1600 

-4500 

-2300 

-2300 

600 

-1200 

-1500 

-2900 

-1200 

-1200 

900 

- 600 

- 600 

-4500 

-2300 

-2400 

1300 

-1300 

-1600 
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Table A-i , Comparison of Viking 1 and Viking 2 Measured VLCA 

Pore os to Preflight Analytical Predictions (Continued) 


Member /Moa s . Mo . 


Minimum 

( Compression), lb, 


Maximum 
(Tension), lb. 


Predicted Viking 1 Viking 2 Predicted Viking f Vi.king .. 2 _ 
Stage I Burnout/Sta ge 11 Igni tion (Model VllI ) 


750/CY 186S 
751/CY 187S 
752/CY 188S 
753/CY 189S 
754/CY 190S 
755/CY 19 IS 


750/CY 186S 
751/CY 187S 
752/CY 188S 
753/CY 189S 
754/CY 190S 
755/CY 1 9 IS 


750/CY 186S 
751/CY 187S 
752/CY-188S 
753/CY 189S 
754/CY 190S 
755/CY 1 9 1 S 


750/CY 186S 
751/CY 187S 
752/CY 188S 
753/CY 189S 
754/CY 190S 
755/CY 1 9 1 S 


-2600 

-2500 

2000 

300 

-1300 

-1300 

1800 

100 

-2700 

-2500 

2100 

300 

-2500 

-2500 

1900 

400 

-1400 

-1300 

2000 

300 

-2700 

-2600 

2000 

600 

Stage II Burnout (Model 

JV) 


-1400 

-1400 

400 

0 

- 700 

- 700 

400 

- 100 

-1400 

-1400 

400 

- 100 

-1300 

-1500 

400 

0 

- 700 

- 700 

900 

100 

-1400 

-1600 

900 

100 

Centaur MES 

II (Model 

I) 


- 700 

- 900 

0 

0 

- 500 

- 500 

100 

0 

- 900 

- 900 

0 - 

0 

- 800 

-1000 

0 

100 

- 500 

- 500 

0 

0 

-_80.0 

-1000 

0 

100 

Centaur MECO 

II (Model 

JV) 


-1600 

-1500 

1300 

400 

i - 800 

- 800 

900 

400 

i -1600 

-1600 

1400 

600 

1 -1500 

-1600 

1600 

800 

) - 800 

- 900 

900 

300 

) -1700 

-1900 

900 

700 


Note: Compression = Negative (-) > Tension .^—Positive (+) 

All values in the above. table have been rounded off to ±100 lbs. 
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APPENDIX B 


Comparison of Measured and An. .. i- . ' v edicted 
Member Forces for the Vova craft 


Estimated flight loads for the Voyager spacecraft were obtained by applying 
the measured Launch Vehicle/Voyager interface accelerations to the base of the 
spacecraft analytical model. The resulting loads are comparable to those 
obtained from the launch vehicle/payload composite model under the 
corresponding environments. 

Since the interface accelerations play a very prominent role in the shock 
spectra approach, the flight measurements will be compared with those obtained 
analytically from the transient analysis. Figure B-l shows this comparison 
for the launch event. The analytical interface acceleration was obtained by 
applying the forcing function representing the lateral overpressure condition 
to the Launch Vehicle/Payload composite model. The amplitude of the 
analytical interface acceleration is somewhat higher than those of 
flight-measured values. However, due to different frequency content, it is 
not certain that the analytical values will produce higher loads. Similarly, 
Figure B-2 shows the interface acceleration comparisons for the Stage I 
Burnout (STG I BO) event. Here not only the amplitudes of the analytical and 
flight data are similar but also the frequency contents are characteristically 
very close. It should be noted that some of the forcing functions used in the 
transient analysis are synthesized based on the experiences from the previous 
flights, and certain conservatism has been built into these forcing 
functions. Vet the resulting analytical interface acceleration is not much 
greater than that of the flight measurements. 

Next, the f Tight loads will be compared with the corresponding design loads. 
Since only the Stage 0 Ignition (launch) and STG I BO events were considered 
in the loads analysis, the comparison will be made by first listing the shock 
spectra design loads for the two events, then listing the ratio of flight 
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loads to the corresponding design loaos as shown in Table B-l. Most of the 
design loads are more than twice the flight loads (the flight to design load 
ratio is less than 0.5). This confirms the design postulations that the shock 
spectra approach will provide conservative loads and in view of the 
uncertainties, the conservatism is reasonable. However, it must be emphasized 
that the comparisons are made for the launch and STG I BO events only. 

Similar comparisons are made for the transient analysis predicted loads as 
shown in Table B-2. In this comparison, the ratios of flight to predicted 
loads are much larger than the ones in the previous comparison. In fact, some 
members have the ratio greater than 1.0 which means that the flight loads are 
greater than the corresponding transient analysis predicted loads. One may 
observe that the transient load analysis indeed does provide more accurate 
loads prediction than tha-t of the shock spectra approach, of course at a 
higher cost. It should also be noted here that the transient leads analysis 




TIME, s 

figure B-2. Interface’ aecelerat Lons for Stage I burnout 

Figure B-2. 

for Voyager was performed not for the purpose of obtaining the design loads 
but rather for the verification of the shock spectra loads. Had the transient 
loads analysis been used for design purposes, a loads analysis factor would 

have been used to multiply the resulting loads to provide more conservative 
design loads. 




Table R-l. Comparison of Maximum Flight: and 
Shook Spectra Design Loads 


Member 


71917 
c « 71927 

£ £ 71937 

4 | 71947 

•H <U 

* £ 71957 

< 71967 


Design 


(lb) 


Launch 


Fliglit/Deslgn 


Stage I Burnout 

Design 

Flight/Design 

(lb) 

A 

B 

820.0 

0. 19 

0.26 

770.0 

0.24 

0.29 

850.0 

0.31 

0.23 

990.0 

0.23 

0. 19 

810.0 

0.25 

0.22 

630.0 

0.27 

0.28 


1000.0 

980.0 
1020.0 
1020.0 

730.0 

710.0 



80017 

80016 

80106 

80107 

80177 

80176 


620.0 

760.0 

460.0 

900.0 

740.0 

530.0 


30007 

30017 

30027 

30047 

30057 

30067 

30077 

30087 


1750.0 

1850.0 

570.0 

1910.0 

1900.0 

1310.0 

1320.0 

450.0 


2400.0 

2540.0 

140.0 

2650.0 

2600.0 

980.0 
1050.0 

200.0 




















































Solid Rocket High-Gain 

RTG Latch Truss Motor Ring Antenna Truss 


labu> 15-;.. Comparison of Maximum Flight and Transient 
Prod io led Loads 


LAUNCH 


Predicted Predicted 


MAX a. 


71917 

71927 

71937 

71947 

71957 

71967 


347.6 

413.6 
563.2 

565.5 

492.6 
313.4 


493.7 

192.6 

257.4 

454.4 
343.0 

466.7 


30007 1021.0 
30017 940.3 
30027 56.1 
30047 1321.0 
30057 1005.0 
30067 592.2 
30077 770.0 
30087 124.9 


0.35 0.58 

0.26 0.76 

0.29 0.63 

0.40 0.79 

0.40 0.64 

0.41 0.70 


0.54 0.78 

0.47 1.68 

0.48 0.67 

0.72 0.94 

0.94 1.19 

0.35 0.78 


0.54 1.10 

0.61 0.96 
1.18 2.16 
0.49 0.88 

0.54 0.87 

0.20 0.76 

0.32 0.79 

0.42 0.99 


344 . 0 

316.5 
360.9 
396.3 

316.6 
269.5 


450.6 

164.2 

350.1 

344.4 
315.9 

441.4 


694.2 

673.3 
54.8 

928.2 

890.9 

553.7 

564.5 

106.4 


Flight . 

Predicted 

A 

B 

0.43 

0.70 

0.33 

0.68 

0.33 

0.61 

0.40 

0.72 

0.43 

0.74 

0.57 

0.81 

0. 66 

0.81 

0.53 

0.79 

0.22 

0.40 

0.98 

1.12 

1.08 

1.21 

0.24 

0.26 

0.80 

1.09 

0.89 

1.17 

1.47 

2.50 

0.63 

0.89 

0.71 

0.91 

0.26 

0.61 

0.30 

0.70 

0.68 

1.42 


STACK 1 BURNOUT 


Flight 

Predicted Predicted 


212.3 
306.7 

423.3 
450.9 
354.5 
144.0 


489.2 
169.9 
176.7 

610.2 
582.2 
324.0 


0.73 1.00 

0.60 0.73 

0.62 0.46 

0.51 0.38 

0.57 0.50 

1.18 1.23 


0.98 0.85 

0.76 1.08 

0.96 1.04 

0.94 0.86 

0.94 0.88 

0.51 0.34 


897 , 

,4 

0 , 

.96 

0 . 

.91 

925 , 

.5 

0 . 

.96 

0 . 

,93 

22 , 

.6 

1 . 

,74 

2 . 

,11 

1080 . 

.0 

0 . 

.91 

0 . 

,86 

1088 . 

,3 

0 . 

.88 

0 . 

86 

506 . 

,3 

0 . 

.48 

0 . 

76 

470 . 

,3 

0 . 

,67 

0 . 

56 



























Scan Platform 

Latch Truss Science Boom 


Table' 15-2. (continued) 


LAUNCH 


STACK i BURNOUT 


Flight 

Predicted P r e d i c t e d 


FI: eht 


Member 


48027 

48037 

48047 

48057 

48067 

48077 

48087 


66.3 

240.8 

268.9 

341.4 

448.5 
187.3 
201.0 
178.2 


338.7 

487.1 
488.9 

259.1 


0.50 1.13 

0.49 0.88 

0.53 0.79 

0.61 0.89 

0.59 0.92 

0.44 1.08 


291. 

7 

183. 

2 

514. 

8 

342 

3 

734. 

.9 

253, 

.8 

366, 

. 1 


0.57 0.96 

0.67 1.10 

0.71 1.11 

0.62 1.13 

0.53 1.13 

0.62 0.95 

0.57 1.19 


59.0 

196.4 

220.6 

255.1 

361.2 
97.2 

104.3 

128.4 


215.4 
153.8 
403.6 

233.5 

587.5 

253.5 
305.0 


286.5 

408.2 
407.0 

233.2 


Predi 

cted 

A 

B 

0.57 

0.92 

0.57 

0.77 

0.62 

0.82 

0.33 

1.00 

0.69 

0.93 

0.64 

0.93 

0.64 

1.10 

0.48 

0.85 

0.82 

0.97 

0.87 

0.96 

1.07 

1.12 

0.97 

1.12 

0.75 

0.92 

0.72 

0.77 

0.70 

0.90 

1.00 

1.08 

0.72 

1.01 

0.77 

0.95 

1.14 

1.24 


Plight 
Predicted Predicted 


212.6 

230.2 

305.2 

373.4 

505.0 

394.4 

423.1 
434.0 


0.30 0.30 

0.82 0.75 

0.76 0.71 

0.91 0.85 

0.88 0.85 

0.41 0.35 

0.41 0.35 

0.35 0.41 


0. 

41 

0. 

45 

0. 

00 

0. 

98 

0. 

.95 

0. 

OO 

OO 

0. 

,90 

0. 

,84 

0. 

CO 

OO 

0. 

,82 

0, 

.99 

0, 

.99 

0, 

.43 

0, 

.39 
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Tnblo 15-2. (continued) 

MAX <t q 

Ight Flight 

let od Predicted Predicted 


STACK I BURNOUT 
Flight 

Prod ie ted Prod ic t ed 


68011 

3234.0 

68021 

3753.0 

68031 

2507.0 

68041 

2143.0 

68051 

3398.0 

68061 

3774.0 

68071 

2056.0 

68081 

2038.0 


0.38 

0.30 

0.26 

0.29 

0.34 

0.40 

0.27 

0.37 


0.73 

0.75 

0.51 

0.80 

0.80 

0.96 

0.71 

0.60 


2632.0 

3248.0 

2922.0 

2798.0 

3087.0 

3188.0 

2944.0 

2335.0 


0 . 

29 

0 . 

27 

0 . 

30 

0 . 

26 

0 . 

27 

0 , 

,24 

0 , 

.23 

0 

.35 


0.48 0.28 

0.34 0.24' 
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APPENDIX C 


Comparison of Acceleration Levels Obtained 
Using the Generalized Shock Spectra Method 
To Those Obtained From the Mass Acceleration Curve 


Figure C-l shows the comparison of acceleration levels obtained from the mass 
acceleration curve to those obtained from the generalized shock spectra method 

for two spacecraft, Galileo (GLL) and the International Solar Polar Mission 
(ISPM) spacecraft. 


DYNAMIC RESPONSE 
ACCELERATION, G 



Figure C-l. Comparison of acceleration levels 


•! 

i 

3 

1 




- 87 - 


APPENDIX D 


A Sample Problem for the Recovered 
Transient Analysis 


The proposed payload transient analysis technique will be demonstrated on a 
realistic complex structural system, namely, the International Solar Polar 
Mission (ISPM) spacecraft. The ISPM consists of two separate spacecraft, one 
sponsored by NASA and the other by the European Space Agency (ESA). Figure 
D-l shows the schematic of the ISPM spacecraft with the major structural 
components identified. Table D-l contains the descriptions of the 
mathematical model. 



F i ^ure D- l . 


ISPM spacer raft 


Table D-l. 
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J.SPM Finite Element Model 


Number of Grid Points 

Number of start,. Degrees of Freedom 

Number of Finite Elements 

Number of Dynamic Degrees of Freedom 

Number of Vibration Modes of Interest 

Frequency ^ 8etal “ d “ Load « Ankysis 


NASA S/C 

112 ’ 
400 
298 
150 
51 


ESA S/C 


122 

430 

221 

87 

84 


Total 

234 

830 

519 

237 

135 

40 

143 Hz 


U ; i;. “ r.u rm r 

buttle launch vehicle configuration consist of e ^ SP * Ce 

and hydrogen tank (ET) and two tw obiter, external oxygen, 

y gen tank (tT), and two solid rocket boosters (SRBs). For 

i nterp 1 anetary trajectories an tik u ■> i a . - 

payloads to reach " 

composite system in the lift-off c-fiouration tw , “WShuttle/IUS 
Which the Shuttle structural system is designed J .ZrollTclZ T °" 

t„rrt: b,ter 

descent, TAEM (terminal arp- SOn * Space °P era tion, entry and 

*« «« JrcTrrs 

:: b ;r:: f ::: r«r n9 events - ° f " -“itr that 

n the 1 f t-off environment will be considered. 

Prior to the ISPM project, another planetary spacecraft, Galileo (GUI has 
been designed to be launched by the identical Shuttle/IUS launch Lte 
system. The Shutt 1 o/IUS/Ga li leo composite model has been n e 1 

dynamic environments representing the lift-off and ah * • * * f th ® 

GLL model. The results of th ! and 3bort using a simplified 

IV.B.l and IV B 2 for the 1 3,13 ^ 1 ^ US6d aS described in Section 

th. analysis of the Shutt le/IUS/ ISPM system. The first 
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Figure D-2. STS/ IUS with ISPM payload 


composite system, Shu it 1 e/ 1 US/ Galileo, is -represented by 150 normal modes at 
^he interface DOF and their corresponding frequencies, i.e., and 
sj* The modal response time histories |q(t)| were also available. The 
cantilevered ISPM model is represented by 40 normal modes and their 


v 


w* 
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correspond ing frequencies, i.e., and j\|J, A constant modal 

damping cycle c/c = 0.01 is assumed for all the modes. First, a mini 
loads analysis, in which tne ISPM model will be subjected to a base motion 
obtained from the IUS/Galileo interface accelerations, will be performed. The 
mini loads analysis is similar to the rigid body interface acceleration method 
of Section IV.B.l without the correction term for the elasticity of the 
spacecraft. The results will be compared with the recovered transient 
analysis. The comparisons will be made on the member loads of the selected 
components. 

Figure 0-3 shows the I'JS/Galileo interface acceleration for the lift-off 
condition. These accelerations have been used as forcing functions for the 
mini loads analysis. Since the ISPM spacecraft is in the same class as that 
of the Galileo spacecraft, the approximate dynamic environment for the ISPM 
can be estimated from Figure D-3. In this case, approximately a three G peak 
acceleration will be experienced at the base of the ISPM spacecraft in the 
longitudinal direction and two G peak acceleration in the lateral direction. 
Table 0-2 shows the loads obtained by the recovered transient method for the 
NASA/ESA spacecraft adapter truss. The comparison of these loads to those 
obtained by the mini loads analysis shows that reasonable agreement has been 
achieved. Since the spacecraft adapter truss consists of structural members 
connecting the two spacecraft, the loads are mainly derived from the low 
frequency motion of the interface. Therefore, it is not surprising that the 
mini loads analysis predicts the loads accurately. However, for the 
structural members dominated by the high frequency local motion, the 
discrepancies between the results of these two analyses can be substantial. 
Figure D-4 shows the time histories of the loads obtained by a mini loads 
analysis and recovered transient analysis for a RCS thruster outrigger. These 
loads are order of magnitude different from each other. 
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